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Abstract. This contribution describes a statistical model for decaying quantum 
systems (e.g. photo-dissociation or -ionization). It takes the interference between 
direct and indirect decay processes cxplicitcly into account. The resulting expressions 
for the partial decay amplitudes and the corresponding cross sections may be considered 
a many-channel many-resonance generalization of Fano's original work on resonance 
lincshapes [Phys. Rev 124, 1866 (1961)]. 

A statistical (random matrix) model is then introduced. It allows to describe 
chaotic scattering systems with tunable couplings to the decay channels. We focus 
on the autocorrelation function of the total (photo) cross section, and we find that it 
depends on the same combination of parameters, as the Fano-parameter distribution. 
These combinations are statistical variants of the one-channel Fano parameter. It 
is thus possible to study Fano interference (i.e. the interference between direct and 
indirect decay paths) on the basis of the autocorrelation function, and thereby in the 
regime of overlapping resonances. It allows us, to study the Fano interference in the 
limit of strongly overlapping resonances, where we find a persisting effect on the level 
of the weak localization correction. 
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1. Introduction 

Molecular photodissociation [HEIIHISEIIE] an d atomic autoionization [Zj are examples of 
quantum-mechanical decay processes: Within the dipole approximation, the absorption 
of a photon excites the quantum system into an energy region, which allows the 
dissociation and/or ionization of the system. Technically speaking, one observes the 
decay of an initial wave packet \a) within a scattering system, via accessible open 
channels. Given a scattering system with Hamiltonian H, the partial decay amplitudes 
t a (E) are of central importance, as they define the partial and total (photo) cross 
sections: 

t a (E) = (¥ a \E)\a) a a (E) = a \t a (E)\ 2 a(E) = J> c (£) . (1) 

c 

Here, ty( a \E) is the post- controlled scattering solution of the problem. This means 
that, asymptotically, there is one single outgoing plane wave in channel a. The sum 
in the definition of the total cross section o~(E) runs over all open channels. The 
proportionality constant a"o depends on the excitation mechanism. As we will not discuss 
those mechanisms, we set o~o = 1. This description reduces the whole process to a half 
collision processes [HI El CDj. Note that we will mostly use the terminology from photo 
absorption, for simplicity, being well aware of the fact that other excitation mechanisms 
are equally possible. 

We aim at a description of chaotic scattering systems, which have a large interaction 
region (in terms of elementary Planck cells), but a finite number of open channels. Such 
a situation occurs frequently in the photodissociation of simple molecules, where the 
potential surface allows trajectories coming from infinity to enter the interaction region 
only through so called "bottlenecks" (transition states). In this situation, scattering 
trajectories may have very short dwell times: They either do not find the bottleneck to 
enter the interaction region, or, due to symmetry, leave the interaction region after only 
a few bounces (Ehrenfest time). These are direct processes. The remaining trajectories 
"randomize" and leave the interaction region only after the average dwell time, which 
is proportional to the size of the openings. They contribute to the indirect decay. 

In [TT] IT2*] . Alhassid and Fyodorov formulate and solve a statistical model for 
the total cross section. This quantity has the convenient property that it can be 
expressed in terms of the Greens function Hlj . With the help of Feshbach's 
projection formalism f5J Ei], the full Hilbert space is divided into an orthogonal 
sum of a number of quasi bound states, and an M-channel continuum. In an 
analytical expression for the autocorrelation function of o~(E) is derived, by using the 
supersymmetry method. The result is quite similar to the "Verbaarschot-Weidenmiiller- 
Zirnbauer" (VWZ) integral [Hj. 

In particular dynamical systems, deviations from the statistical model are quite 
possible. In the semiclassical regime, for instance, Policot-Ruelle resonances may give 
rise to characteristic signatures in the autocorrelation function of the total photo cross 
section [THJ [TH] . A different source for deviations, is the interference with direct decay 
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processes (Fano interference). In the limit of isolated resonances, it leads to asymmetric 
lineshapes |2T)ll21| . The statistical model of Alhassid and Fyodorov has been generalized 
recently, to include direct processes [221 EHj The basic formalism has been developed 
in nuclear reaction theory [21] (and references therein); see also reference |2*o] . 

Fano interference (Beutler-Fano lineshapes) has also been studied in mesoscopic 
transport. It is shown in reference j2H] that Fano interference can be observed in the 
conductance of quantum dots coupled to single mode leads. Even though these results 
are not simply related to the present model, they still point at important applications. 
Quite generally, a Fano resonance may be used as a kind of built-in interferometer, 
which, for instance, is very sensitive to decoherence effects |27j . Also temperature ef- 
fects as well as Coulomb interaction effects can be studied |28j . 

In this paper, we address the question how to detect and quantify Fano interference 
in the many channel case, and in the regime of overlapping resonances. For that 
purpose, we review in section 121 the statistical model for half collision processes, as 
developed in |2~2"l l2"3*j . The mapping formalism of reference [22], together with an 
important generalization is presented in section El We then consider the limit of 
isolated resonances, where the resonance lineshapes can be described by a complex 
Fano parameter (section EJ). The statistics of the Fano parameter has been considered 
in [2H], for the one channel case. Here, we identify those quantities, which determine 
the Fano-parameter distribution in the many channel case. For each channel, one can 
find a statistical analog of the one-channel Fano parameter. In sectional we turn to the 
study of the autocorrelation function. Surprisingly, it depends on very similar quantities. 
This allows to quantify Fano interference, beyond the regime of isolated resonances. In 
section El we study Fano interference in the "semiclassical regime" , where many channels 
are strongly coupled to the interaction region. Even there, we find a noticeable effect. Its 
strength is comparable to the weak localization correction in mesoscopic transport [30J. 
A summary is given in section [7| 

In the present paper, the discussion is restricted to the case of orthogonal invariance, 
which allows to use the VWZ-integral to compute autocorrelation functions. For 
broken time-reversal symmetry, the corresponding result for the S-matrix element 
correlations has been given recently [HI]. It would allow to repeat the analysis of total 
cross section autocorrelation functions for that case. 

2. The model 

We follow the approach adopted in [T2j [221 [2Hj where resonant and direct processes are 
treated separately by means of the Feshbach approach [T^l EE]: the Hilbert space of 
excited states is divided into a subspace of bound states (Q-space), and its complement. 
The Schrodinger equation restricted to that complement is treated as an auxiliary 
scattering system. Choosing the bound space appropriately, it looses all resonant 
features. Its scattering solutions serve as a basis, in which the resonant scattering 
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via the subspace of bound states is formulated. 

We consider the scattering problem (or the half collision problem) in a small energy 
interval, away from any thresholds, such that the number of open channels, denoted by 
M, remains constant. Then, the partial decay amplitudes, defined in equation (0), may 
be written as 

t a (E) = (a|« out ) + (a| Ft 1 ( - in V |a out ) ) H cS = H Q -mVV\ (2) 

& — -Hcff 

where the vector \a) has been decomposed \a) = \a m ) © |a out ) according to the 
partitioning of the Hilbert space, above. Here, Hq is a iV x iV-matrix of the original 
Hamiltonian projected onto Q-space. The matrix V contains the transition amplitudes 
between the basis states in Q-space, and the scattering states |c). It has M column 
vectors (the channel vectors), each of dimension N. 

Note that t a (E) may be considered as a transition amplitude between the state \a) 
and the scattering states \a) of the auxiliary resonance-free scattering system. These are 
typically not the free scattering solutions. To obtain the transition amplitudes into the 
free scattering states, one also needs the scattering matrix of the auxiliary scattering 
system. Below, we focus on the total photo cross section, where such considerations are 
of minor importance. 

To obtain the total photo cross section, one may sum over all partial cross sections, 
or use a particular form of the optical theorem ^3]- It relies on the fact that the 
scattering solutions \ty( c \E)) form a complete basis. Therefore: 



M 

a(E) = V (a|* (c) ) {¥ c) (E)\a) = (a\ 5(E - H) \a) = - Im (a| G{E~) \a) 



c=l 

Ha^f-TT- 1 Im 



( |a in ) + ivr V |a out ) ) f * ( |a in ) - m V |a out ) ) 

E — ti e B 



(3) 



Generally, we will assume that the column vectors of V are pairwise orthogonal. If that 
is not the case, one may use a singular value decomposition [22]: V = Vqs, where s 
is a unitary matrix. It allows to replace V by Vq, whose column vectors are pairwise 
orthogonal. Due to V \a ont ) = Vq |a° ut ) with |a° ut ) = s |a out ), the expression Q remains 
unchanged except for a redefinition of the direct decay amplitudes. 



2.1. The statistical model 

For the derivation of the above expressions for t a (E) and c(E), we have used an arbitrary 
but fixed basis in Q-space. In that basis, we assume V, a m and a out to be fixed, and take 
H from the Gaussian orthogonal ensemble (GOE). For the full collision problem, the 
statistical properties of the resulting S-matrix have been studied in detail f7[ ESI E3]- 
Indirect photo decay has been studied in |T2], and the effects of direct decay paths has 
been considered in j22 122] • 

For the description of decay processes, we have 3M + 1 independent parameters, 
which define the model: 



Statistical quantum decay 



5 



• The norm of the column vectors of V measure the coupling strengths to the decay 
channels. 

• The direct decay amplitudes a™' = (c|a out ) give the overlaps of the initial state 
with the scattering states of the auxiliary scattering system. They describe direct 
decay processes, which may may lead to Fano interference. 

• The overlaps of a m with the column vectors of V, i.e. (c| |a m ), may give rise to 
Fano interference, also. As we will see, this type of interference is not equivalent to 
the one before. 

• The norm of the component of |a in ), which is orthogonal to all column vectors of 
V. That component is responsible for the indirect decay processes. 



3. The mapping formalism 

The mapping formalism has been introduced in [221 • ^ allows to express the partial 
decay amplitudes as well as the total photo cross section in terms of an extended 
scattering matrix. This is convenient if analytical results for the statistical scattering 
model are to be transfered to the present case. Here, we will use it to compute the 
average (this section) and the autocorrelation function (section EJ) of the total photo 
cross section. In its original formulation, |a m ) is required to be orthogonal to the 
column vectors of V. Here, the formalism is generalized, to allow for non-zero overlaps 
between |a in ) and the channel vectors. This leads to the distinction between external 
and internal Fano interference. 

Assume, the initial state |a m ) has overlap with the channel vectors, and that channel 
vectors are coupled to the continuum with non-vanishing transmission coefficients. In 
that case, we can always find a |a 10 ), such that 

\a ln ) = |o>) + vr V \a [o ) , (4) 

where |ag n ) is orthogonal to all column vectors of V. Then, one can write for the 
transition amplitudes t a [equation (j2J)] and the total cross section [equation (j3J]: 

t a = ( a \a out ) + (a\ V ] 1 ( \o$) - m V |a° ut ) ) |c4 ut ) = |a out ) ± i |a io ) 

E — H cS 

a(E) = \\a out \\ 2 - tt" 1 Im ( |o>) + iirV |a° ut ) ) f ( |o>) - m V |a° ut ) ) 

E — H cS 

Following [2Hj, this suggests to introduce the following extended scattering matrix: 
S(E, 6) = 1 1+M - 2m H/t 1 W W = ( o> 5/(2tt) , V ) , (6) 

— fteflf 

which contains only the orthogonal part of |a m ). Now, one has to find transformation 
matrices u(S) and v(S), such that S'(E, 5) = u(6) S(E, 5) v T (5), may be used to construct 
the partial decay amplitudes, and the total photo cross section. One obtains the 



(5) 



Statistical quantum decay 
following answer: 



u(S) 





(a° ut )t 5/2 \ 




( i 


{a°^) T 8/2 \ 


















v(5) = 








1a/ 









) 




I o 





(7) 



The transformation of S(E, 8) is reminiscent of an "Engelbrecht-Weidenmiiller" 
transformation jSH], although for 5 > 0, u(5) and v(S) are not unitary. The partial 
decay amplitudes and the total photo cross section can be obtained from: 



(a\n ont ) 

t a (E)= { -^—l + hm5~ 1 SUE,S) 

Z o— >U 

a(E) = ||a ou i 2 - ^ Re(a° ut |a° ut ) + lim ^ Re[ 1 + ^ (E, 5) ] . 



(8) 
(9) 



The average total cross section As a first application, we will compute the average 
total cross section (o~(E)) in the center of the spectrum. To this end, we use the fact 
that the average extended S-matrix is given by [T7j : 

(S ab (0,5)) = 5 ab \^ A a = 7r 2 p K|| 2 A = ^ po IKH 2 , (10) 



1 + A a " ru 11 u " u 4 ru 11 u 11 ' 

where for a = 1, . . . , M , v a are the respective column vectors of the coupling matrix V. 
The level density is denoted by po, i.e. the average level spacing is A = 1/(N p ). As 
the average S-matrix is diagonal, we obtain: 

\2 



1 — A h 1 



(11) 



where S is the average of the original scattering matrix (of dimension M). From this it 
follows with equation Q: 

(a) = \\a out \\ 2 + po K n f + - Re (a° ut | (5 - 1) |a° ut ) . (12) 



Provided, the components of |a out ) and \a m ) are real, we may proceed a bit further. In 
equation (J3J), |a 10 ) is chosen such that 

(c|a in ) = 7T (v c \a [o ) = ir\\v c \\ (c\a io ) = 



Therefore: 



I „ out II 2 i _ 1 1 _ in 1 1 2 
\a || + Po ||«0 II 



| „ out 1 1 2 , _ 1 1 _ in 1 1 2 
\a || + po ||Q! || 



E 

c 

E 



out 1 2 



. out 1 2 



(c|a 10 ) 



^ i«i n r 



+ Po |<f) 



Po /A c (c|a in ) . (13) 
A, 



1 + A C 

A c 
1 + A C 



(14) 



On the level of the average total cross section, it makes no difference, whether a m has 
overlap with the external channel region, or internal states (inside Q-space) which are 
connected to the channel region via transmission coefficients. For later use, we define 



the direct photo cross section as cr di] 



a 



out || 2 



and the indirect photo cross section as 



CTind = PO IK 1 ! 2 - 
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4. Fano lineshapes: the limit of isolated resonances 

The effect of direct reaction paths is most evident in the limit of weak coupling, where 
individual resonances can be observed. In this regime, the resonances have asymmetric 
lineshapes (Beutler-Fano profiles) [201 El], i.e. the cross section near resonance can be 
parametrised as 

<T[E) = <7 dir = (Tdir 5— : e = /r> • (15) 



S 



1 mr e 2 + l Yj/2 



While in the one channel case, the Fano parameter q may be assumed real, one needs a 
complex Fano parameter q = qi + i<?2 in the case of many channels. The imaginary part 
of the Fano parameter lifts the resonance, such that the cross section at its minimum 
is no longer zero. As shown below, the background cross section (the cross section at 
large |e|) is <7di r = ||« out || 2 , independent of the imaginary part of the Fano parameter. 

In the limit of isolated resonances, the effective Hamiltonian H e g can be assumed 
diagonal in the eigenbasis of H . Using that basis in equation (J3J), we obtain: 

(a, + i bj)* {a,j — ibj) 



a(E) = ||a out || 2 - - Yl Im 



7T 

3 



E — Ej + iTj/2 



(16) 



where aj = (j|« in ), bj = n J2 c Vj C a° nt and F,- = 27r^ c ^A 2 ,. If these parameters may be 
assumed real (e.g. due to time reversal invariance), we obtain in the vicinity of the j'th 
resonance: 

a(E) « aj(E) = \\a^\\ 2 - ^ Im[ (e - i) (a, - ib f } . (17) 

This expression can be put into the form of equation (fi3|) with a^r = ||a out || 2 . This 
yields the following values for the real and the imaginary part of the Fano parameter: 



2*4 ^.fi + ^Wi-l 



where Dj = it ||a out || 2 Yj/2 > 6 2 due to the Schwarz inequality (it makes sure that q2 is 
always real). 

In the case of a regular cross section, it makes sense to study the individual 
resonances and their lineshapes. However, if the lineshape (i.e. the Fano parameter) 
fluctuates randomly from resonance to resonance, a statistical analysis is more 
appropriate. Let us rewrite equation (|TH|) in terms of M + 1 independent normalized 
Gaussian random variables. We will take into account that |a m ) may have some overlap 
with channel vectors (the column vectors of V) . 



V jc -> S c X c dj = (j\a n ) + 71 Y V jc «c° -* S \X + 7T ^ 



so 



s c a° r nt „ 2 ||« out|12 



b,=nJ2^^ ^ = - 2 ^^E^ s c- (19) 

c u c 

Besides a^h = ||a out || 2 and cr in< i oc Sg, the following parameter combinations can be 
identified: s c , s c a° ut /s , and s c a c ° / s , where c = 1, . . . , M. As these parameters come 
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along with centered Gaussian random variables, their sign is irrelevant. We thus prefer 
to consider their squares: 

2 _ \\ a II 2 _ /y2v 

6 - at S c — \ V jc) 



N c x JC/ 4tt 2 iVpo 



s 2 |a° ut | 2 |a° ut | 2 T c _ r c out s 2 \a™\ 2 \a™\ 2 T c 



— (20) 
s 2 4 7T 2 a ind 4tt 2 s 2 4 vr 2 a ind 4vr 2 ' 

As we will see below, the autocorrelation function depends almost on the same set 
of parameters. That means, an analysis of the autocorrelation function gives much the 
same information, with the advantage that its application is not restricted to the regime 
of isolated resonances. 



One channel case In that case b 2 = Dj so that the Fano parameter becomes real: 

a ~ — - W) a - (21) 

where we have again separated that part of |ct m ) which is orthogonal to V and that 
which is parallel. For precisely that case, the distribution of qi has been computed 
in [2H]. The result is a shifted Lorentzian with width g determined by: 

2 (CtKT) 4 q ind J_ 

- 2 |a out | 2 (V. 2 ) |a out | 2 T x r? nt ' 



7T- 



whereas the shift is given by 



„,OUt T-OUt 



(23) 



That part of the work in j2H|, which deals with time reversal invariant systems, is hence 
fully contained in the present model. 



5. Autocorrelation function of the total photo cross section 

In section 01 direct decay processes have been studied, as they affect the lineshapes of 
isolated resonances. Once the resonances start to overlap, such an analysis is no longer 
possible. Here, we study the interference effects on the autocorrelation function. This 
quantity has the advantage, that there are no restrictions on the coupling strengths to 
decay channels (transmission coefficients). This allows to study Fano interference even 
in the regime of strongly overlapping resonances (section IHJ). 



Internal versus external Fano interference Here, we will study the autocorrelation 
function of the total photo cross section (j3J), with the help of the mapping formalism 
and the VWZ-integral ^7j. The latter provides an analytical expression for the 
correlation function between two matrix elements of the (extended) scattering matrix in 
equation ©. With this, it is possible to compute correlation functions between partial 
decay amplitudes, as well as between total cross sections. It is also possible to choose 
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different initial states a, if that would be of interest. In the case of two total cross 
sections, we would consider the correlation function 

C[a 1} a 2 ](w) = (a 1 (E-w A/2) a 2 (E + w A/2)} -(a^E)} (a 2 (E)} , (24) 

where (. . .) denote a spectral ane/or ensemble average, as appropriate, and A the mean 
level (resonance) spacing. Here, (u^) denotes the total photo cross section for the 
initial state \a\) {\a 2 )). Usually, we find it more convenient to consider the correlation 
function in the time domain: 

C[a ll a 2 ]{t) = j dw e 2niwt C[a u a 2 ](w) . (25) 

Due to the fact that C[o"i, <J 2 ](w) is defined on the unfolded energy axis, t measures time 
in units of the Heisenberg time. 

For simplicity, and because this case is typically considered in the literature, the 
following discussion will be restricted to the autocorrelation function of the total photo 
cross section. Using for a the expression Q from the mapping formalism, (section 02) 
one obtains: 

C[a] = lim I C[S' 00 + S£, S' 00 + S'* ] = 2 lim 1 Re C[S' 00 , S'* ] , (26) 

and for its Fourier transform (see reference [HE]): 

C[a]=Km±C[S' m ,S'* ]. (27) 

The auxiliary scattering matrix S'(E, 5) is defined within the mapping formalism [see 
equations © and Q]. For its matrix element S' 0Q , one obtains: 



M , 2 M 

S'oo = Soo + g 5 J2 t («-'): + («?)« S ^ ] + J E S °* ■ ( 28 ) 

a=l a, 6=1 

This gives C[cr](t) in terms of a linear combination of correlation functions between 
S-matrix elements, where the VWZ-integral |T7] applies. Due to equation (|A.1|) many 
such correlation functions vanish. After a bit of algebra, one arrives at the following 
result: 

6[<r] = oZ d l\ 2A -| EaUr a A a +4n 00 + - £ [r a + r* a + r" + r a + ] II 



2 

a=l 



Oa 



4 E [^T 6 + + r a r 6 ] n afe l (29) 



iW 
a, 6=1 

r a = («° ut ) a (af); T a /a- md r± = | (af% \ 2 T a /a ind . 
If the Fano interference is purely external, then r a = : real, and we obtain the result 



from [2*3] : 



{2 M 1 M 

(2A -^a=l T aK) +4n 00 + 2^r a n„ fl + - Yl T ^ b Tlab 
a=l a, 6=1 



(30) 
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where r c = |a° ut | 2 T c /o~ ind = r° ut as given in equation (|20|). The parameters r c and rf 
are the statistical variants of the Fano parameter; see section 0J Eventually, we call 
them "statistical Fano parameters" for short. If Fano interference is purely internal, 
then r± = K | 2 T a /V ind = r^°, r a = -t±. Therefore: 

(2 A + i Eili ^ A a ) + 4 n 00 + - £ Il ab . (31) 

a,6=l J 

6. Limit of many open channels 

We consider two cases. The first is the absorptive limit. There, all transmission 
coefficients go to zero, while the number of channels goes to infinity: M — > 00. The 
limits are taken in such a way that the sum over all transmission coefficients remains 
finite: ^ C T C = T sum . In that case, the resonance widths stop fluctuating jHZj. In the 
second case, we assume that all transmission coefficients are equal to one, while the 
number of channels is large but finite. We refer to that the semiclassical regime. 

Absorptive limit In this case, the autocorrelation function can be computed in closed 
form (all integrals can be solved, see |Appendix A ), and one obtains: 



tfM = °Le-' J -' J> o r[l-& 2 (i)] + 2 + 0)j , (32) 

where the constant Cq depends on the statistical Fano parameters. Hence, in the absence 
of internal Fano interference, one can obtain a perfect exponential autocorrelation 
function by properly tuning ^ a r a — > 4. 

Semiclassical regime For simplicity, we focus on the two extreme cases. In the first 
case, the decay is dominated by indirect processes, such that the autocorrelation function 
becomes: 

C , M(t)=4 4 d XW(A 2 + n o), (33) 
while in the second case, it is dominated by direct processes. This results in 

C[a] (t) = al d Ea ' fe(T ^ + + 7aTb) 1^ n n . (34) 

As can be seen from equation ([A. 1)1 . the autocorrelation function is then proportional to 
the autocorrelation function of a diagonal S-matrix element (note that all transmission 
coefficients are equal to one). Thus, very large statistical Fano parameters (i.e. the 
dominance of the direct processes) result in a photo cross section which has the same 
statistical properties as a full scattering cross section er to t(.£7) [HH]. Below, we will thus 
compare C[a] (t) for indirect decay with Cfcrtot] (t) for direct decay. Note that for large M, 
one expects the autocorrelation function in both cases to be dominated by an exponential 
decay with e~ Mt 
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Figure 1. The autocorrelation function Cfcr]^) versus Mt, in the absence of direct 
processes (indirect cross section) compared to the same quantity where the direct 
processes dominate. The curves are computed numerically, on the basis of the VWZ 
integral, see | Appendix A| All transmission coefficients are taken equal to one, and the 
number of channels is varied between M = 16 (solid line), 8 (dashed line), 4 (dotted 
line), and 2 (dash-dotted line). The thin dotted line shows the exponential decay: 
2 exp(-Mt). 



In figure Q we simply compare the resulting autocorrelation functions, varying the 
number of channels. A semilog plot is used, and on the ordinate Mt is given, such 
that differences to the exponential decay are easier recognized. The autocorrelation 
functions are normalized, such that C[cr](0) = C , [cr tot ](0) = 2. In the case of indirect 
decay, the autocorrelation function lies above the purely exponential decay (thin dotted 
line), approaching the exponential as M increases. In the case of direct decay, the 
autocorrelation function approaches the exponential from below. 

In figure El the same autocorrelation functions are shown, again. However, here, we 
divide by the classical expectation exp(— Mt) and we plot the autocorrelation functions 
versus time t (in units of the Heisenberg time). In this way, differences to the classical 
expectation are strongly enhanced. For indirect decay, the autocorrelation function has 
been computed semi- classically in [U], with the result: C[a](t) = 2 (1 + 1) e~ Mt . For 
the autocorrelation function of total scattering cross sections, a similar result is not 
available (to the best of the author's knowledge). One may however find a connection 
to the weak localization correction of mesoscopic transport jHOj, because the integral 
over the autocorrelation function in the time domain, gives the average value of the 
modulus squared of the corresponding S-matrix element [0*0] : 

/ dt C[a](t) = (\S aa \ 2 ) = . (35) 
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t 



Figure 2. The same autocorrelation functions as in figure^ but divided by e~ Mt , and 
plotted versus t. Again, the number of channels is varied between M = 2 (dash-dotted 
lines), M = 4 (dotted lines), M = 8 (dashed lines), and M = 16 (solid lines). The 
autocorrelation functions for indirect decay have a positive slope at t = 0, while in the 
opposite case, they have a negative slope. 



M 


a 


b 


c 


2 


1.99934(54) 


0.076(48) 


0.1(15) 


4 


2.00018(35) 


1.991(30) 


3.67(87) 


8 


2.00029(44) 


5.987(37) 


5.0(10) 


16 


2.00057(21) 


13.953(17) 


8.74(49) 



Table 1. The coefficients a,b,c in equation Ij36(l as obtained from a polynomial fit to 
the C[a}(t) as plotted in figure The digits in brackets give the estimated error from 
the x 2 -fit [SI- The first column gives the number of channels. 



To study the behavior of the autocorrelation functions C[<r](t) (indirect decay) and 
C'[c r tot](i) more quantitatively, we perform a polynomial fit to the numerical results, as 
shown in figure 121 To this end we use the following general expression: 

C[a](t) ~ (2 + at + bt 2 + ct 3 + . . .) e~ Mt . (36) 

The results of the fits are given in the tables Q and El We used a polynomial of fifth 
order for the fit, but the values for the higher order coefficients had too large errors to 
be of any use. On the basis of the values found, one may conjecture: 

C[a](t) ~ (2 + 2t+ (M-2)t 2 ) e~ Mt (37) 
C[a tot ](t) ~ (2 - At + (M + A)t 2 ) e~ Mt . (38) 

In the case of indirect decay, we obtain agreement with ^H] up to the linear term. 
Higher order semiclassical corrections to the indirect photo cross sections have not been 
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M 


a 


b 


c 


2 


-4.00109(70) 


10.101(62) 


-27.1(19) 


4 


-4.00039(48) 


12.028(41) 


-32.5(12) 


8 


-3.99930(30) 


15.945(25) 


-46.21(72) 


16 


-3.99951(21) 


23.948(18) 


-77.58(49) 



Table 2. The coefficients a, b, c in equation as obtained from a polynomial fit to 
C[atot](t) as plotted in figure [5] The digits in brackets give the estimated error from 
the x 2 -fit E2- The first column gives the number of channels. 

considered, so far. In the second case, we may actually integrate equation (|3~%j) to recover 
the weak localization correction: 

p*<W«>-4-s* + S5 + " (39) 

in agreement with an expansion of equation (|33|) in inverse powers of M. This may of 
course just be a coincidence. For a definite answer, one should try to integrate X^ M ^ iln 
in closed form. 

7. Conclusions 

This work has started out from a recent formulation of a statistical model for half 
collision processes [221 I2H|- The principal aim was to investigate the possibility to 
detect and quantify Fano interference in cases where the total photo cross section shows 
irregular behavior (that may be randomly fluctuating lineshapes and/or overlapping 
resonances). The analysis of the Fano-parameter distribution j2Hl I2H] m &y be 
appropriate, as long as the resonances are well isolated. However, the cross section 
autocorrelation function, can always be analysed. Surprisingly, it provides practically 
the same information (statistical Fano parameters) as the analysis of the Fano-parameter 
distribution. This allowed us to investigate Fano interference in the regime of strongly 
overlapping resonances. There, Fano interference effects are much weaker (on the level 
of the weak localization correction). 

One might have doubts, whether Fano interference can also be detected in practice, 
when many channels are strongly coupled to the interaction region. Given a dynamical 
system, the main question is, how much averaging over different initial conditions and/or 
different samples is necessary in order to obtain the desired information. To clarify that 
point, numerical simulations with the open kicked rotor [2111111312] are currently under 
investigation. 

Fano interference has turned into a versatile tool in mesoscopic physics, where it 
is used to measure dephasing and decoherence times [2H1I2Z1- So far, only single Fano 
resonances in single mode systems have been applied. It would be desirable, to be 
able to measure dephasing and decoherence times also in multi-mode systems, from 
irregular cross sections. As a first step in that direction, one should include dephasing 
and decoherence into the statistical model. 
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Appendix A. The Verbaarschot-Weidenmuller-Zirnbauer integral 

In the case that Hq is taken from the GOE, the correlation function between two 
scattering matrix elements from equation (jHJ) is given by the VWZ-integral [H| . Fourier- 
transformed into the time domain, it reads 



C[S ab ,S* cd ](t) =l{45 ab 6 cd T a T c A a A c + 2{6 ac 5 bd + 6 ad 6 bc )T a T b n ab } . (A.l) 
Here, the T a are the transmission coefficients, and we used the following abbreviations: 

(t-r)(r+l-t) l-T e (t-r) 



I = dr du 



dr I 

ix(0,t-l) Ji 

r + T a x t — r \ 2 2r + l 



(2u + l)(t 2 - r 2 + x) 2 A = l + 2T e r + T 2 



i \/l + 2T P r + Tfx 



A a = Jl-T a - ■ — " +- — * = ^ 



l + 2T a r + T 2 x l-T a (t-r)J 2u + 1 

n T a T fc x 2 +[T a T b r + (T a + T 6 )(r + l)-l]x + (2r + l)r (t-r)(r + l-t) 



1 + 2T a r + T 2 x)(l + 2T b r + T 2 x) [1 - T a (t - r)] [1 - T 6 (t - r)] ' 

Absorptive limit There, all transmission coefficients go to zero, while the number 
of channels goes to infinity: M — > oo. Both limits are taken in such a way that 
Tq = T sum remains finite. In that case, one find: 

A a -> A = t 

Tl ab -> n 00 = (t - r)(r + 1 - t) - x + (2r + 1) r (A.3) 



r 



1 _ e - T -'2<«0 *°> = / dr / du " ^ + 1 " 



= j( o) noo = i (A 4) 



'max(0,t-l) JO (2« + 1) (t 2 - r 2 + X) 2 

The integrals J (0) and J (0) n 00 can be calculated in closed form This gives: 

1 j(o) n - 1 

It 2- J n °°-2 
For the GOE spectrum, the two point form factor b 2 (t) is given by 

l-b 2 {t) = 2t-t \n(2t + 1) + 9{t - 1) [2-2t + t ln(2t-l)] . (A.5) 

This gives for the correlation function: 

C[S abl S*J(t) = e~ Tst [5 ab 5 cd T a T c (1 - b 2 (t) ) + (5 ac 6 bd + 5 ad 5 bc ) T a T b ] . (A.6) 

The semiclassical regime There, all transmission coefficients are set equal to one 
Vc : T c — 1. Typically the number of channels M is assumed to be large but finite. 
However, the following relations really hold for arbitrary M. In that case, one finds: 

T a ,T b ^0 



A = t :T a ^0 
Al =0 :T a ^l ab 




T a -*0,T b ^l ,(A.7) 
T ,T 6 - 1 
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where n o is given in equation (jA.3|) . and 

^ ^ x + r t — r , . . 

Hoi =t U u = — - - + — 7 • A.8 

1 + zr + x 1 + r — t 

Finally, 

1 -. Z<"> = f dr f d„ («-)(r+l-«) / 1+r-t \ » 

The resulting integrals are computed numerically. 
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